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We examine the effects of an unparticle IA as a possible source of missing energy in the 
p-wave decays of B meson. The dependence of the differential branching ratio on the 
Kq (Ki) —meson's energy is discussed in the presence of scalar and vector unparticle op- 
erators and significant deviation from the standard model value is found after addition of 
these operators. Finally, we have shown the dependence of branching ratio for the above 
said decays on the parameters of the unparticle stuff like the effective couplings, cutoff scale 
An and the scale dimensions du- 



I. INTRODUCTION 



Flavor changing neutral current (FCNC) processes induced by b — > s transitions are not allowed 
at tree level in the Standard Model (SM), but are generated at loop level and are further suppressed 
by the CKM factors. Therefore, these decays are very sensitive to the physics beyond the SM via 
the influence of new particles in the loop. Though the branching ratios of FCNC decays are 
small in the SM, quite interesting results are obtained from t le experiments both for the inclusive 
B -> X s £+£- Q and exclusive decay modes B -» K£+l~ 41 and B -► K*£ + i~ 
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These 
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results are in good agreement with the theoretical estimates 

Among different semileptonic decays induced by b — ► s transitions, b — ► svv decays are of 
particular interest, because of the absence of a photonic penguin contribution and hadronic long 
distance effects these have much smaller theoretical uncertainties. But experimentally, it is too 
difficult to measure the inclusive decay modes B — > X s vv as one has to sum on all the X s 's. 
Therefore, the exclusive B — * K (K*) vv decays play a peculiar role both from the experimental 
and theoretical point of view. The theoretical estimates of the branching ratio of these decays are 
Br (B — * Kvv ~ 1(T 5 ) and Br(B -» K* vv ~ 10 6 ) [al whereas, the experimental bounds given 

n n 

by the -B-factories, BELLE and BaBar, on these decays are [101. Illl]: 



Br (B -> Kvv) < 1.4 x 10~ 5 
Br (B -» K*vv) < 1.4 x 10~ 4 . 



(1) 



These processes, based on b — > svv, are very sensitive to the new physics and have been studied 



extensively in the literature in the context of large extra dimension model and Z' models fl^ . [13I . 
Any new physics model which can provide a relatively light new source of missing energy (which is 
attributed to the neutrinos in the SM) can potentially enhance the observed rates of B — > K (K*) + 
missing energy. Recently, H. Georgi has proposed one such model of Unparticles, which is one of 
the tantalizing issues these days [141 ]. The main idea of Georgi's model is that at a very high energy 
our theory contains the fields of the standard model and the fields of a theory with a nontrivial 
IR fixed point, which he called BZ (Banks-Zaks) fields [ijj. The interaction among the two sets is 
through the exchange of particles with a large mass scale My. The coupling between the SM fields 
and BZ fields are nonrenormalizable below this scale and are suppressed by the powers of My. The 
renormalizable couplings of the BZ fields then produce dimensional transmutation and the scale 
invariant unparticle emerged below an energy scale Ay. In the effective theory below the scale Ay 
the BZ operators matched onto unparticle operators, and the renormalizable interaction matched 
onto a new set of interactions between standard model and unparticle fields. The outcome of this 
model is the collection of unparticle stuff with scale dimension dy, which is just like a non- integral 
number of invisible massless particles, whose production might be detectable in missing energy and 
momentum distributions [lfj. 
This idea promoted a 



at colliders 
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ot of interest in unparticle physics and its signatures have been discussed 
17|, Q, Q,|2flin low energy physics [2l|, Lepton Flavor Violation 



221 ]. unparticle 



physics effects in B s mixing [23], and also in cosmology and astrophysics 24j. Aliev et al. have 
studied B — > K (K*) + missing energy in unparticle physics [25] in which they have studied the 
effects of an unparticle Wasa possible source of a missing energy in these decays. They have found 
the dependence of the differential branching ratio on the K (i*T*)-meson's energy in the presence 
of scalar and vector unparticle operators and then using the upper bounds on these decays, they 
put stringent constraints on the parameters of the unparticle stuff. 

The studies are even more complete if similar studies for the p-wave decays of B meson such 
as B -> K% (1430) + ft (ft is missing energy) and B -> K x (1270) + ft, where K% (1430) and 
K\ (1270) are the pseudoscalar and axial vector mesons respectively, carried out. In this paper, we 
have studied these p-wave decays of B meson in unparticle physics using the frame work of Aliev 
et al. 25J. We have considered the decay B —* Kq (K\) vv in SM although for these modes no 
signals have been observed so for, but in future B-factories where enough data is expected, these 
decays will be observed. These SuperB factories will be measuring these processes by analyzing 
the spectra of the final state hadron. In doing this measurement a cut for high momentum on the 
hadron is imposed, in order to suppress the background. Therefore, unparticle would give us a 
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unique distribution for the high energy hadron in the final state, such that in future B-factories 
one will be able to distinguish the presence of unparticle by observing the spectrum of final state 
hadrons in B -> (K, K*, Kq, K x ) + $ [25]. 



The work is organized as follows: In section II after giving the expression for the effective 
Hamiltonian for the decay b — > svv, we define the scalar and vector unparticle physics operators 
for b — > slA. Then using these expressions we calculate the various contributions the decay rates 
of B — > Kq (K\) + fi both from the SM and unparticle theory in Section III. Recently, Grinstein 
et al. gave comments on the unparticle [26f] mentioning that Mack's unitarity constraint lower 
bounds on CFT operator dimensions, e.g dy > 3 for primary, gauge invariant, vector unparticle 
operators. To account for this they have corrected the results in the literature, and modified the 
propagator of vector and tensor unparticles. We will also give the expressions of decay rate using 
these modified vector operators in the same section. Finally, section IV contains our numerical 
results and conclusions. 

II. EFFECTIVE HAMILTONIAN IN SM AND UNPARTICLE OPERATORS 

The flavor changing neutral current b — * svv are of particular interest both from theoretical and 
experimental view. One of the main reason of interest is the absence of long distance contribution 
related to the four-quark operators in the effective Hamiltonian. In this respect, the transition to 
neutrino represents a clean process even in comparison with the b — * S7 decay, where long-distance 

n 

contributions, though small, are expected to present [2?]. In Standard Model these processes are 
governed by the effective Hamiltonian 



H eff = ^^V tb V*C 10 sr (1 - 75) hv ltl (1 - 75) v (2) 

where VtbV t * are the elements of the Cabbibo-Kobayashi Maskawa Matrix and C\q is obtained from 
the Z° penguin and box diagrams where the dominant contribution corresponds to a top quark 
intermediate state and it is 

_ D{x t ) . . 

sin o w 

8 W is the Weinberg angle and D (xt) is the usual Inami-Lim function, given as 
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with xt = rnl/rriy^. 

The unparticle transition at the quark level can be described by b — > sU, where one can consider 
the following operators: 



Scalar unparticle operator 



HA 1V U 

Vector unparticle operator 



Cs^-s-y^Ou + C P ^s W bd»O u (5) 



C Vj^^ b °u + C A^Sl^O» u . (6) 



HA 2 HA 
The propagator for the scalar unparticle field can be written as 0,0 

A d 



17] 



J d'xe^ (0 \T0 U (x) O u (0)| 0) = ij—^ (-P 2 )^ 2 (7) 



with 



lOvr 5 / 2 T(d u + 1/2) 
d » {27T ) 2d uT(d u -l)T(2d u Y U 

III. DIFFERENTIAL DECAY WIDTHS 

In Standard Model the decay B — > Kq (K\) + E is described by the decay B —* Kq (K\) vv. 
At quark level this process is governed by the effective Hamiltonian defined in Eq. (J2J) which when 
sandwiched between B and Kq {K\) involves the hadronic matrix elements for the exclusive decay 
B — > Kq (Ki) vv. These can be parameterized by the form factors and the non- vanishing matrix 

H 

(Kl (pf) \s W b\ B(p)) = -i [/+ (q 2 ) (p +p') ii + /_ {q 2 ) g„] . (9) 

where = (p + p') . Using the above definition and taking into account the three species of 
neutrinos in the Standard Model, the differential decay width as a function of Kq energy (Ek*) 
can be written as [271 ]: 

rTTSM r 2 n 2 / : 

v„ - " '"IS 



elements for B — > Kq are 



with A (a, b, c) = a 2 + b 2 + c 2 - 2ab - 2bc - 2ca and q 2 = M| + M£. - 2M B E K * . Here /+ (q 2 ) 
and /_ (q 2 ) are the form factors which are the non-perturbative quantities and can be calculated 
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using some models. The one we have used here was calculated by using Light Front Quark Model 
(LFQR) by Cheng et al. [27( and these can be parameterized as: 

F(a= m 5 

V ' 1 - o? 2 /Ml + b (? 2 /M|) 2 
and the fitted parameters are given in Table HI 

TABLE I: The parameters for B — > Kq form factors. 





F(0) 


a b 


/+ 


-0.26 


1.36 0.86 


f- 


0.21 


1.26 0.93 



Similarly, for B — > K\ transition the matrix elements can be parametrized as 28] 



(K 1 (k,e)\V^\B(p)) = is^Ms + M Kl )V x {q 2 ) 

V 2 (q 2 



(K 1 (k,e)\AJB(p)) 



-(p + fc) M (e*-g) 



^'M B + M Kl 
2M Kl 



-e* u p a k (3 A(q 2 ) 



(11) 
(12) 



M B + M Kl 

where = s^^b and = 57^756 are the vector and axial vector currents respectively and £* is 
the polarization vector for the final state axial_vector meson. In this case we have used the form 
factors that were calculated by Paracha et al. 



281 ] and the corresponding expressions are: 



A (a) 
Vi(s) 

V 2 (s) 



A(0) 



(1 - s/M B ) (1 - s/M'i) 
l-s/M 2 )(l-s/M' 2 



M 2 - Mf Ci 



(13) 



V2(0) 



2M Kl 



Vb(0) 



1 - s/M B * s/MS. ) M b ~ M Kl (1 - s/Ml) (1 - s/M'i) 



with 



A(0) 
V 2 {0) 

The differential decay rate can be calculated as 

dT SM GW 



•(0.52 ±0.05) 
-(0.24 ±0.02) 
(0.39 ±0.03). 



25[]: 



dEK x 



2 9 7r 5 M| 



V tb V t *f Idol 2 \ l/2 \M SM \ 2 



(14) 



(15) 
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where 
\M SM \ 



2 8g 2 AU(g 2 )| 2 



+ 



{M B + M Kl f M 2 Ki 



V 2 (q 2 



(M B + M Kl y 



+ (M B + M Kl ) 2 (A + l2M 2 Ki q 2 ) |Ux (g 2 



-A (M 2 - - g 2 ) Re(V* (q 2 ) V 2 (q 2 ) + V* 2 (q 2 ) Vi (q 2 )) 



(16) 



and A = A (M|, M 2 K ^q 2 ) with q 2 = M B + M 2 Kx - 2M B E Kl . 

Now in decay mode -B — ► if,* (ifi) + E, the missing energy shown by E can also be attributed 
to the unparticle and hence the unparticle can also contribute to these decay modes. Therefore, 
the signature of two decay modes B — > if<* (ifi) vv and -B — ► if<* (K\)U is required like the one 
done for B -> if (if*) ^ and S -> if (if*) W in the literature [25(] . 



A. The Scalar Unparticle Operator 

Using the scalar unparticle operator defined in Eq. ([5]) the matrix element for B — ► ifgW can 
be written as 



^ = ^<^0OI^m(Cs + CpT&)&|£(p)>^ 



Now the decay rate for B —* ifg^ can be evaluated to be: 

dT su 1 



El*. -Ml. \M SU \ 2 



where 



\M 



SUI 2 



1 Pl A 2d 



M 2 B + M 2 K *-2M B E K ^" ' 



x /+ (q 2 ) (M 2 - M 2 ;,) + /_ (q 2 ) (m 2 + Af|j - 2M B E K * 
Following the same lines, the corresponding matrix element B — > ifiZ^ is 

^ = a^7 l s "> ^ 5 + Cp75 ) 5 I B (P)> 

= JLc s (e*.g) (Mb + MkJU!^ 2 ) 

-(M b - M^)U 2 (g 2 ) - 2M Kl (U 3 (<? 2 ) - V {q 2 )) 



O u , 



and the differential decay rate is 
dT su _ M B A du 



^ = ^Pt\ C s\ 2 \ V ° (? 2 ) I' ( E k ~ <) 3/2 ( M b + M 2 Kl - 2M B E Kl 



2 N3/2 



(17) 



(18) 



(19) 



(20) 



(21) 
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One can see from Eq. (|18|) and Eq. (|2ip that the scalar unparticle contribution to the decay 
rate depends on Cp, Cs, du and An, therefore one can see the behavior of decay rates for the said 
decays on these parameters which will be hoped to get constraint once we have experimental data 
on these decays. This we will do in a separate section. 



B. The Vector Unparticle Operator 



The matrix element for B — > KqIA using the vector unparticle operator defined in Eq. and 
the definition of form factors given in Eq. ([9]) can be calculated as: 



M V K H = -^ I (KZ(p')\s ll ,(C v + C Al5 )b\B(p))O i 



1 



J^CaIU (q 2 ) (p +P% + /- (q 2 ) qM. 



The differential decay rate is then 
dT vu 1 A. 



dEx* 



8ir 2 m B A 



Jfc* \C A \ 2 \f+ {q 2 ) | 2 [m% + Ml* - 2M B E K * 



du-2 



' E h~ M k 



X < 



M 2 B + M 2 K *+2M B E K *^ + 



(m% - Ml* 



M% + M 2 - 2M B E K * 



For B — > K\ case the matrix element for B — > K\U is 
M ™ = I ^T(Ki{p')\~ s lACv+CA^)b\B(p))O i 



A§^iK (M B + M Kl ) Vl (q 2 ) - i (p + p% (e* • q) 



2M^, , 2 , 2 „ C A 2A(g 2 ) 



(■^af3£ P P 



(22) 



(23) 



(24) 



and the differential decay rate will be: 
dT vu 1 A th , 



dEx 1 8TT 2 mB A 



M 2 



(M B + M Kl )V 

(M B + M Xl ) 4 (3M^ + 2M 2 M 2 ri - §M B M 2 Ki E Kl + M 2 ^) |^ (g 2 )| 2 

+2M B - M 2 Kx ) | V 2 (q 2 ) | 2 + 4 (M B + M Kl ) 2 
(Mb^ - M|J (M 2 Kl - E 2 Kl ) M 2 B (VxV 2 * + V 2 V 1 * 



(25) 
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The total decay width can be obtained if we integrate on the energy of the final state meson in the 



range M K(Kl) < E K{Kl) < ( Af| + M 2 ) /2M B for B - K {K x ) + E. 



261 ] in which they have men- 



Recently, Grinstein et al. have given comment on the unparticle 
tioned that Mack's unitarity constraint lower bounds on CFT operator dimensions, e.g. du > 3 
for primary, gauge invariant, vector unparticle operators. To account for this they have corrected 
the results in the literature, and modified the propagator of vector and tensor unparticles. The 
modified vector propagator is 



d 4 xe lPx (0 \T (0£ (x) O u (x)) | 0) = A du {-gT + aP»P»/P 2 ) (P 2 ) 



2\<4/-2 



(26) 



Here P is the momentum of the unparticle, Ad u is defined in Eq. ([8]) and o / 1( in contrast to 
the value a = 1 which was considered by Georgi [14]) but is defined as: 



2 {d u - 2) 

a = 



W=W (27) 

By incorporating this factor a in the vector unparticle operator the Eqs. (|23p and (|25[) get modi- 
fication and the modified result of the decay rate for B — > KqU is 



dT 



vu 



dE K * 



l^l 2 \U if) I 2 (Ml + Ml* - 2M B E Kt ) du ~ 2 y/Ek-M*. 



\U(i 2 )\ 



I 



\ 



\ a ( M B- M k 

M B + M 2 K5 +2M B E K5 j + 



M 2 B + M£. - 2M B E K * 



+ |/_ [q 2 ) | 2 (a - 1) [m 2 b + M 2 K , - 2M B E K » 
+2(a-l) {f + {q 2 )f„{ q 2 ))(M 2 B -M 2 K , 



(28) 



Similarly, for B 
dT vu 



K\IA the result becomes 
1 Ada 



dE Kl ' 87r 2 m B A 2d "" 



E 2 Kl ~ M 2 Ki (q 2 ) 



\M U \ 2 + \M 22 \ 2 + |-M 33 | 2 + \M U \ 2 + \M 23 \ 2 + \M 2i \ 2 + \M 34 \ 2 



(29) 
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with 



\M U 



\M 



1 



22 



\C 



V 



M 2 Ki (M B +M Kl ) 2 q 2 
(M B + M Kl f (3M 2 Kl (M 2 B + M 2 Kl - 2M B E Kl ) - a (M 2 B M 2 Kl - M 2 B E 2 Kl )) \Vi (q 2 ) f 



\M 



33 



\Cv\ 



\M 



■l-i 



M 2 Ki (M B +M Kl ) 2 q 2 
M 2 (E 2 Ki - M 2 Kl ) (a (M 2 - M 2 Ki f + (2M B E Kl ) 2 - (M 2 + M 2 K f) \V 2 (q 2 ) \ 

2 1 



\Cv\ 



M 2 Ki (M B +M Kl ) 2 q 2 



\M 



23 



4M 2 (M B + M Kl ) 2 (E 2 Ki - MfJ (a - 1) M 2 Ki \V 3 (q 2 ) - V [q 



\Cv\ 



2M 2 



\M 



21 



(M B +M Kl ) 2 q 2 

M 2 B (M B + M Xl ) 2 (E 2 Kl - M 2 Kl ) (M 2 B + M 2 Kl - 2M B E Kl -a(M 2 B - M 2 Kl )) 
(Vi(q 2 )V 2 *(q 2 )+V 2 (q 2 )V 1 *(q 2 )) 

\Q |2 1 

Vl M 2 Ki (M B +M Kl ) 2 q 2 

2M Kl (M B + M Kl f ((1 - a) M 2 B (E 2 Ki - M 2 Ki )) (V x (V 3 - V )* + (V 3 - V ) V?) 



\M 



3 1 



\Cv\ 



2M Kl (M B + M Kl ) (M 2 B - M 2 Kl ) 



M 2 Ki (M B + M Kl ) 2 q 2 
M 2 B (E 2 Ki - M 2 Ki ) (a - 1) (V 2 (V 3 - VoT + (V3 - V ) V 2 *) 



(30) 



One can easily see that Eqs. (p8|) and (|29|) reduces to the Eqs. (p3|) and (p5|) respectively, if one 
puts o=l. 

IV. RESULTS AND DISCUSSIONS 

In this section we present our numerical study for the B — > (ifi) + ^ where we try to 
distinguish the unparticle physics effects from that of the SM. In Standard Model E which is 
the missing energy is attributed to the neutrinos where as in the case under consideration, this is 
attached to the unparticle. Therefore the total decay rate can be written as 

r = r 5M + r". (31) 
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Here V s M is the Standard Model contribution (B — ► Kq (K-\) vv) where as T u is from the unparticle 
(B -> (ifi) W) to the decay £ K% (K x ) + In ref. [25] it is pointed out that the SM process 
B ^ K (K*) vv provides a unique energy distribution spectrum of final state hadrons and present 
experimental limits on the branching ratio of these processes are about an order of magnitude below 
the respective SM expectation values. They have used experimental upper limit on the branching 
ratio of B — > K (K*) vv decay to estimate the constraints on the unparticle properties. 
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FIG. 1: The differential branching ratio for B — » Kq (K\) + $ as a function of hadronic energy Ek* {Eri) 
is plotted. Top panel is for B — > Kq+ Jfi and bottom is for B — > Ki+ ~$. The other parameters are dy — 1.9, 
A w = 1000 GeV, C7 P = C s = 2 x 10" 3 and C v = C A = 10~ 5 . Solid line is for SM, dashed line is for scalar 
operator and long-dashed line is for the vector operator. 



In case of B — > Kq {K\ ) vv there is no experimental limit on the branching ratio of these 
decays, but these will be expected to be measured at future SuperB factories where they analyze 
the spectra of final state hadron by imposing a cut of on the high momentum of hadron to reduce 
the background. To calculate the numerical value of the branching ratio for B — * Kq (K\) vv in 
SM we have to integrate Eqs. (|10p and (|15p on the energy of the final state hadron. Thus after 
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integration, the values of the branching ratios in SM are: 

Br (B — > KqVv) = 1.12 xl(T 6 (32) 
Br (B -> K\vv) = 1.77 x 1(T 6 

With these values at hand, we have plotted the differential decay with for B — > Kq {K\) + E 
as a function of the energy of the final state hadron Ek* (Eki) an d by fixing the parameters 
of unparticle from ref. [25[] in Fig{T] One can easily see from the figure that the signature of 
unparticle operators are very distinctive from the SM for the final state hadron's energy. Just like 
B — > K (K* ) + E the distribution of unparticle contribution is quite different when we include a 
vector operator (a = 1) for the highly energetic final state hadron. For the other values of a we 
will discuss this issue separately. Thus the Super B-factories will be able to clearly distinguish the 
presence of unparticle by observing the spectrum of final state hadrons in B —* Kq {K\) + E in 
complement to B —> K (K*) + E. 

In Fig. 2 and Fig. 3 we have shown the sensitivity of the branching ratio on the scaling 
dimension du for different values of the cut-off scale Au by using the same values of Cs, Cp, Cy 
and Ca as we have used for Fig. 1. We can see from this figure that the branching ratio is very 
sensitive to the variable du and Ay. The constraints on the vector operator are more stronger then 
the scalar operators and constraints for B — > E are better then the B — * K\-\- E decays. 

After showing the dependence of branching ratio on du and Au what we have shown in Fig. 
4 is the sensitivity of the branching ratio of B — ► ^ W1 th the effective coupling constants 

of scalar and vector unparticle operators. One can see that B — > Kq+ scalar unparticle operator 
shall constrain the parameter Cp and B — ► Kq+ vector unparticle operator shall constrain the 
parameter Ca- Thus observing this decay we can get some useful constraint on Cp and Ca which 
provides us the signature about the unparticle physics. Similarly, we have shown the dependence 
of the branching ratio of B — > K\+ E on the effective coupling constants in Fig. 5. It is shown 
that if we consider the scalar operator then only dependence is on Cs, whereas if we consider the 
vector operators then the decay rate depends both on Cy and Ca- 



As we have a 
the unparticle 



ready mentioned that, Grinstein et al. have recently given their comment on 
261 ] mentioning that one regards Mack's unitarity constraint lower bounds on CFT 
operator dimensions, e.g., du > 3 for primary, gauge invariant, vector unparticle operators. To 
account for this they have corrected the results in the literature, and modified the propagator 
of vector and tensor unparticles. The modified expressions of decay rate for the processes under 
consideration are given in Eq. (|28p and Eq. f|29|) . To incorporate this modification in vector 
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FIG. 2: The branching ratio for B — ► Kq+ Tfi as a function of c?^ for various values of Ajj- Top panel is for 
the scalar operator and bottom is for the vector operator. The values of coupling constants is same as taken 
for Fig. 1. Solid line is for An = 1000 GeV, dashed line is for An = 2000 GeV and long-dashed line is for 
An = 5000 GeV. The horizontal solid line is the SM result. 



unparticle operator, what we have shown in Fig. 6 is the fractional error 



dr \ l dV 



T dE K * , K > / I V dE K * , K j 

A = ^ ) J a=x ^ /a (33) 

dr 



TdE H^) /a=l 

where the difference is between the spectrum of B — » Kq (K\ ) U using vector unparticle operator 
with a = 1 and a = 2 (du — 2) / (dy — 1) with 3 < dy < 3.9. It is clear from the graph that if we 
increase the unparticle scaling dimensions du the contribution of vector unparticle operator to the 
decay rate decreases significantly because of the increase in the inverse powers of cutoff scale An 
(see Eqs. Qgfy and flggj). 

Just to conclude: The study of these p-wave decays of B mesons will not only provide us 
information about SM but it also indicate the physics beyond it and in future, when enough data 



13 



0.0001 



bo 

.5 0.00008 



? 0.00006 
t 0.00004 \ 
W 0.00002 \ 

o \ — 



1.4 



1.5 



1.6 



1.7 



1.8 



1.9 



0.0006 




FIG. 3: The branching ratio for B — > Ki+ ^ as a function of cfo for various values of Ay. Top panel is 
for the scalar operator and bottom is for the vector operator. The values for coupling constants is same as 
taken for Fig. 1. Solid line is for An = 1000 GeV, dashed line is for Au = 2000 GeV and long-dashed line 
is for Au = 5000 GeV. The horizontal solid line is the SM result. 



is expected from the Super B factories, we believe that these decays will take us step forward to 
the study of unparticle as a source of missing energy in flavor physics. 
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FIG. 6: Fractional error A in the spectrum for the decay B — > Kq (K\) + vector unparticle operator as a 
function of energy of final state hadron. Top panel shows the value for B — » Kq and the bottom panel is 
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is for du = 3.2, dashed line is for du = 3.4 and dashed-double dotted is for d u = 3.6. 
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